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Sensing magnetic nanoparticles using nano-confined ferromagnetic 
resonances in a magnonic crystal 
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We demonstrate the use of the magnetic-field-dependence 
of highly spatially confined, GHz-frequency ferromagnetic 
resonances in a ferromagnetic nanostructure for the detec¬ 
tion of adsorbed magnetic nanoparticles. This is achieved 
in a large area magnonic crystal consisting of a thin fer¬ 
romagnetic film containing a periodic array of closely 
spaced, nano-scale anti-dots. Stray fields from nanopar¬ 
ticles within the anti-dots modify resonant dynamic mag¬ 
netisation modes in the surrounding magnonic crystal, 
generating easily measurable resonance peak shifts. The 
shifts are comparable to the resonance linewidths for high 
anti-dot filling fractions with their signs and magnitudes 
dependent upon the modes’ localisations (in agreement 
with micromagnetic simulation results). This is a highly 
encouraging result for the development of frequency- 
based nanoparticle detectors for high speed nano-scale 
biosensing. 

Magnetic biosensors, in which biological analytes are tagged 
with magnetic nanoparticles (MNPs), have excellent potential 
for solid-state point-of-care medical diagnosticsMIl. The tech¬ 
nique is intrinsically matrix-insesntive^, can compete with 
industry-standard immunoassays® and can be combined with 
magnetic separation methods®. The central element of a mag¬ 
netic biosensor is a detector for the stray or ‘fringing’ mag¬ 
netic fields generated by magnetised MNPs which are used 
to label, typically in-vitro, analytes of interest within a bi¬ 
ological sarmle. Previously used sensors include SQuIDs®, 
Hall sensors®, ferromagnetic rings^® and magneto-impedance 
devicesCSl. However one of the most widely used methods 
is tha t em^plo ying magnetoresistive (MR) magnetic field sen- 
sors®®2KHI which are typically fabricated with at least one 
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lateral dimension on the order of 10-100 /im. An MNP is de¬ 
tected when its stray (or ‘fringing’) magnetic field modifies 
the quasi-static magnetic configuration in the ferromagnetic 
MR device. This changes the device’s resistance, enabling 
electronic MNP detection. 

It can however be challenging to minimise noise in conven¬ 
tional MR sensors when reducing the sensor size due to ther¬ 
mal instabilities of the device’s magnetic configurationSSEl 
A suggested approach to overcome this is to use intrinsically 
high frequency detection methods exploiting the strong field 
dependence of resonant magnetisation dynamics®^ which 
can be reliably driven in isolated nanostructures. These dy- 
nam ks c an be driven electrically in spin torque oscillators 
(e g | i8 | 2H l24] therein) which have been predicted to 

retain high field sensing signal to noise ratios at sub-100 nm 
dimensions!Furthermore, real time electrical detection of 
the dynamics^SEll will pave the way for high speedE^l, nano¬ 
scale MNP sensing for solid-state fiow cytometry^^®!] 

In this work we use a large areal^^Sl magnonic crystal 
(MC)IS1 to macroscopically probe resonant GHz-frequency 
magnetisation dynamics which are spatially confined to nano¬ 
metric regions and experimentally demonstrate their use for 
MNP sensing. In contrast to continuous ferromagnetic layers, 
nanostructured MCs (e.g. Fig.[T^) exhibit a number of distinct 
ferromagnetic resonance (FMR) modes with different lateral 
spatial localisations within the crystal’s nano-periodic struc- 
tureEEl. A previous study demonstrated that FMR modes 
within anti-dot-based MCs are sensitive to magnetic nanos¬ 
tructures fabricated within the anti-dots, leading to a pro¬ 
posal for a magnonic biosensor to detect captured MNP-based 
tagsESl. Here we demonstrate that stray magnetic fields gen¬ 
erated by captured MNPs within an anti-dot-based MC do in¬ 
deed generate clear resonance peak shifts which are mode- 
dependent and which approach the resonance linewidth when 
the fraction of MNP-filled anti-dots is high. Note that reso¬ 
nances within MNPs can be be detected directly however only 
broad, relatively weak signals have been observed previously 
for small collections of MNPs®. Rather, in this work, we 
detect MNP-induced changes to very well defined resonances 
within a periodically nanostructured, high quality ferromag¬ 
netic film. 
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Fig. 1 (a) Scanning electron micrograph of a portion of a magnonic 
crystal (MC). The scale bar is 1 long, (b) Schematic of the 
experimental setup showing the MC plaeed faee down above a 
micro-stripline. 


Large area MCs (4x4 mm^) consisting of arrays of 0.3 
jum wide anti-dots with edge-to-edge spacings of 0.15 
were fabricated in a 30 nm thick N i 8 oFe 2 o ferromagnetic layer 
using deep ultraviolet lithographyEEl (Fig. [^). The MCs’ 
FMR modes were probed at room temperature using broad¬ 
band, stripline-based, magnetic field modulated FMR spec¬ 
troscopy (Fig.[^). This technique measures the derivative of 
finite width FMR peaks with respect to the external magnetic 
field, F/ext, at a fixed GHz frequency. MC spectra are measured 
both before and after the addition of cluster-shaped MNPs 
(diameters ^ 0.1 —0.3 ^m). Two micromagnetic simulation 
methods have also been used: (i) a time domain (‘ringdown’) 
simulation using MuMax3l3SI which subjects the MC to a field 
pulse and exploits Fourier analysis on the resulting dynamics 
to extract frequency-resolved information (e.g.®^), and; (ii) an 
eigenmode method (e.g.®) which directly calculates the sys¬ 
tem’s resonant modes and mode profiles for a given F/ext (car¬ 
ried out using FinMag which is the successor to Nmag®!). See 
the supplementary information for additional details. 

An experimental FMR spectrum for a bare MC (ie. no ad¬ 
sorbed MNPs) obtained at 11.5 GHz is shown in Fig. in 
which two high amplitude resonances near the extremes of 
the measured F/ext range can be identified. The frequency 
(/) dependence of the resonance fields, F^res, of these two 
FMR modes compare well with those predicted from time do¬ 
main simulations for the extended (E) and side (S) modes®® 
(Fig.|^). In Fig.[^ we show Fourier transformed data from 
the time domain simulation obtained at ^o^ext = 200 mT, 
showing the E and S modes together with their spatial locali¬ 
sations inside the MC’s unit cell. The E mode is concentrated 
in bands between rows of anti-dots oriented perpendicular to 
F/ext while the S mode is localised between neighboring anti¬ 
dots (see also the schematic in Eig.[^). In both simulation and 
experiment, a number of lower amplitude modes lie between 
the side and extended modes, demonstrating good agreement 
in terms of the overall mode structure however these modes 
will not be discussed in this communication. 

MNPs were adsorbed onto the MCs at F/ext = 0 by plac¬ 
ing 12 juL of diluted MNP suspension on the MC’s surface 
which then was allowed to dry in ambient conditions be- 
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Fig. 2 (a) Experimental, field-resolved FMR trace (/ = 11.5 GHz) 
showing differential absorption peaks corresponding to FMR modes 
in the MC. (b) Comparison of experimental and simulated resonant 
frequencies versus /^oF/ext for the side (S) and extended (E) modes, 
(c) Frequency-resolved, Fourier transformed time domain simulation 
data at A^gF/ext = 200 mT for a single unit cell of the MC with and 
without a 150 nm MNP at the centre of the anti-dot. The Fourier 
amplitude has been differentiated with respect to / to facilitate 
comparison with experimental spectra. Insets show the localisation 
of the resonant dynamics for the S and E modes with lighter shading 
indicating a stronger dynamic magnetisation component 
perpendicular to //ext- The anti-dot boundary is shown as a blue 
circle, (d) Schematic showing the spatial localisation of the E and S 
modes together with a MNP and a sketch of its stray magnetic field. 
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Fig. 3 (a) SEM image of MNPs on the MC with a 300 nm wide 
scale bar. Out of 33 holes counted, 30 contain MNPs. In 7 of these, 
the ‘captured’ MNP extends outside of the anti-dot. There are 5 
isolated nanoparticles on the MC’s upper surface. FMR traces 
obtained at 11.5 GHz showing the side mode (b) and extended mode 
(c) before (‘bare’) and after the addition of MNPs. Multiple traces 
have been taken following repeated removal and replacement of the 
MC, a requirement for MNP application, enabling an estimation of 
the associated experimental uncertainty in 0.5 mT). Traces have 
been vertically scaled and vertically shifted to locally normalize the 
differential absorption signals to that obtained with a bare MC. 


fore re-measurement. Scanning electron microscopy of dried 
MNPs on a MC reveals irregularly shaped MNPs with the ma¬ 
jority lying inside the anti-dots (see Fig. and caption). The 
MNPs generate an upward shift in iTres for the side mode res¬ 
onance (Fig.[^) and a downward shift for the extended mode 
resonance (Fig. [^). Both shifts exceed experimental uncer¬ 
tainty (see caption of Fig. [^. 

The observed shifts can be qualitatively understood by con¬ 
sidering the field generated by an idealized magnetized MNP, 
located at the center of an anti-dot and magnetized along F/gxt 
(Fig. igi). Treating the MNP as a dipole, the y-component 
of its stray magnetic field will, to a first approximation, rein¬ 
force F/ext at E mode region but oppose F/gxt at the S mode re¬ 
gion. Thus, for a given experimental measurement frequency, 
a larger F/gxt must be applied to attain the side mode reso¬ 
nance condition when MNPs are within the anti-dots. Simi¬ 
larly, the extended mode will be observed at a lower external 
field. These predictions are consistent with the experimental 
results (Fig. [^,c) as well as analogous numerical results ob¬ 
tained for anti-ring structures®. 

To verify these arguments more rigorously, simulations 
were repeated with a 150nm wide spherical MNP within the 
anti-dot (see supplementary information for further details in¬ 
cluding the F/exfdependent MNP moment). The time domain 
simulation result for the MC in the presence of a MNP with 
I^qHqxi = 200 mT is shown as a red dotted line in Fig.[^. The 
S mode’s frequency is indeed decreased, consistent with that 


part of the MC experiencing a lower net field. Likewise, an 
increased resonance frequency is predicted for the extended 
mode. To extract numerical values for these shifts, eigenmode 
simulations were carried out at both /ioF/ext = 37 mT (to probe 
the side mode as per Fig J%) and at 206 mT (to probe the ex¬ 
tended mode as per Fig.|^). This yielded MNP-induced fre¬ 
quency shifts of -f0.311 GHz for the E mode and -0.085 GHz 
for the S mode, in good agreement with time domain simu¬ 
lations. The local gradient of the /res vs. F/es data (Eig. |^) 
was then used to convert the frequency shifts into equivalent 
field shifts. This yielded —6.9 mT for the extended mode (45 
GHz/T) and +2.0 mT for the side mode (43 GHz/T). Although 
the observed shift will depend on the MNP coverage, these 
simulations which assume one 150 nm wide MNP per anti¬ 
dot correctly predict the order of magnitude and sign of the 
resonance field shift (Pig. where, albeit, not every hole is 
filled and there is a distribution of MNP sizes). 

To study MNP coverage effects, we measured a second MC 
with an equivalent anti-dot lattice geometry. We carried out 
consecutive applications of diluted solutions of MNPs with in¬ 
creasing concentration, c, obtaining PMR traces and imaging 
the MC via SEM before and after the application of each solu¬ 
tion. Representative SEM images are shown in Pigs. & -d for 
each concentration. Increases in c visibly increase the MNP 
coverage, resulting in an increased peak shift for both the ex¬ 
tended and side modes. Por the lowest coverage, a shift of the 
extended mode is measurable and just higher than the experi¬ 
mental uncertainty (Pig.|^). At the highest coverage however 
the resonance peak shift is much larger and approaches the 
peak-to-peak resonance linewidth. There was some variation 
in the measured shifts for different / however no clear, widely 
applicable trends could be determined (except at low / for the 
side mode, and thus low F/ext^ where a smaller shift presum¬ 
ably resulted from a lower F/ext-induced MNP moment). In 
Pigure|^, we have averaged the S and E modes’ peak shifts 
over the measured frequency range (11.5 — 16 GHz) and plot¬ 
ted the averaged shifts versus c. Reliably measurable shifts 
of the E mode are observed over the entire frequency range. 
However, a higher concentration must be used to register a 
consistent shift for the side mode. Notably, the simulated 
shifts of +2 mT for the S mode and —6.9 mT for the E mode 
compare well with the shifts observed at c = 0.225 jug/juL 
where we are close to having every anti-dot filled with a MNP 
(Pig. 1^)) and thus closest to the simulation condition. The 
MNPs also generate a coverage-dependent linewidth increase 
with a ^1.5 times increase on average in the peak-to-peak res¬ 
onance linewidths at the highest c (Pig. |^). However, the 
linewidth broadening does not dominate the observed shifts 
in that the minimum of the differential absorption line consis¬ 
tently moves in the direction of the peak shift (Pig.[^). This 
is qualitatively consistent with that expected for a collection 
of differently sized MNPs all acting in unison with the shift 










magnitude depending upon the size of the MNP. 
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Fig. 4 (a-d) 2.25 wide images of the MC following consecutive 
applications of aqueous MNP solutions of increasing concentrations, 
c (jjg/jiiL). (e) Differential absorption peak for the extended mode at 
12 GHz in the bare MC and following adsorption of MNPs for 
increasing values of c. Extended and side mode peak shifts (f) and 
relative peak-to-peak linewidth increase (g) as a function of MNP 
solution concentration. Error bars combine both the uncertainty in 
the shift measurement at each frequency and the spread of peak 
shifts over the measured frequency range (11.5-16 GHz). 


Conclusions 

We have used an anti-dot based magnonic crystal (MC) to 
experimentally demonstrate the ability to detect magnetic 
nanoparticles (MNPs) via their influence on spatially con- 
flned, high frequency, ferromagnetic resonance modes. MNPs 
are preferentially captured within the holes which, depend¬ 
ing upon a mode’s spatial localisation, leads to an increased 
or decreased resonant frequency. Resonance shifts are re¬ 
produced well by micromagnetic simulations and observable 
even for quite low anti-dot flllings (^ 15%). A non-dominant 
linewidth broadening is observed at high MNP coverages. Our 
results are directly applicable to conflned modes in isolated 


spintronicES or magnonic®! nanostructures. This is encour¬ 
aging for the development of frequency-based spintronic de¬ 
vices such as spin torque oscillators for nano-scale magnetic 
biosensingESin applications such as flow cytometry®®!. No¬ 
tably our observed frequency shifts are signifi cantly larger 
than measured spin torque oscillator linewidths®IEM3 how¬ 
ever electrical detection will rely on close proximity of the 
MNP to the sensing layer and high GHz/T field sensitivities. 
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Supplementary information 

Magnonic crystals and magnetic nanoparticles 

The MCs were composed of 30 thick nm Ni 8 oEe 2 o layers 
covered by an Au capping layer (8 nm thick for the data in 
Manuscript Eigs. 2 and 3 and 10 nm thick for the data in 
Manuscript Eig. 4). They were fabricated on Si substrates 
using deep ultraviolet lithography with de positio n via elec¬ 
tron beam evaporation followed by lift-orf®®. The anti¬ 
dot lattice geometries were measured using scanning elec¬ 
tron microscopy (SEM) with antidot diameters and array pitch 
rounded to the nearest 10 nm. MNPs were nanomag-D iron- 
oxide nanoparticles (micromod Partikeltechnologie GmbH) 
with a colloidally stabilized dextran surface, a quoted MNP 
width of 130 nm and a quoted solids content of 25 mg/mL. 
The latter was used to calculate the concentrations of the di¬ 
luted MNP solutions (diluted using e-pure water). Magne- 
tometry on freeze dried MNPs was carried out at 300K using 




























a MPMS3 SQuID magnetometer (Quantum Design Inc.) in 
VSM mode. The magnetic moment per unit volume (Fig.[^ 
was calculated assuming an iron-oxide density of 5.24 g/cnr. 
SEM was carried out with a Zeiss 1555 VP-FESEM and a FEI 
Verios 460 SEM. 



Fig. 5 [Supplementary figure] Magnetic moment per unit volume 
for the MNPs at 300K. 


Ferromagnetic resonance spectroscopy 

FMR modes were probed at room temperature using broad¬ 
band, stripline-based, magnetic field modulated FMR spec¬ 
troscopy (Manuscript Fig. lb). During measurement, the 
MC’s (01) axis was closely aligned with stripline which it¬ 
self was aligned with F/ext^ the latter modulated at 220 Hz. 
The technique uses a interferometric receivei®! and lock-in 
amplifier (SR830, Stanford Research Systems) to measure the 
external magnetic field (F/ext-)derivative of finite width fer¬ 
romagnetic resonances in the sample at a set frequency, /, 
and stepped F/gxt- To obtain the FMR traces, F/gxt was in¬ 
creased in steps over a range typically spanning ^ 0 — 250 
mT with F/ext measured at each step using a FH54 Teslame- 
ter (Magnet-Physik Dr. Steingroever GmbH). The resultant 
FMR spectra were measured both before and after the addi¬ 
tion of cluster-shaped MNPs to the MC’s surface. The MCs 
are characterised by a non-zero remanent magnetisation with 
the stepped F/ext sweeps for the measurements carried out for 
a single field polarity. Thus, the bias field and the spatially av¬ 
eraged y-component of the MC’s magnetic moment remained 
aligned during the FMR experiments. A microscope cover- 
slip between the MC and the stripline was used in all mea¬ 
surements to avoid MNPs, when present, rubbing off onto the 
stripline. A PVC block placed on the stripline board with¬ 
out contact to the stripline itself was used to ensure consistent 
placement of the sample with respect to the stripline. This 
enabled excellent reproducibility even with repeated removals 
and replacements of the MC. Reproducibility was confirmed 


for each measurement with the associated uncertainty in the 
peak position typically being on the order of 0.5 mT. This can 
be seen in Manuscript Figs. 3b,c where we show traces ob¬ 
tained after repeated removals and replacements of the sam¬ 
ple. Upon adding MNPs to the MC, we consistently observed 
an overall decrease in the signal amplitude which was stronger 
than the MNP-induced linewidth broadening effects. This sig¬ 
nal amplitude reduction increased with MNP concentration, 
suggesting a broadband absorption of the microwave power 
by the MNPs. To consistently compare data obtained with and 
without MNPs, we vertically scaled the FMR data so that all 
traces had the same peak-to-peak amplitude. It was also some¬ 
times necessary to introduce a small vertical offset, typically 
on the order of a few tens of /nY at most. 

Micromagnetic simulations 

Two micromagnetic simulation methods were used in this 
work. Both methods simulated a single unit cell of the MC 
with periodic boundary conditions in the x and y directions, an 
11 X 11 Tiled’ macro-geometry® for determining the demag¬ 
netizing field and the following magnetic parameters: damp¬ 
ing a = 0.008, nil intrinsic anisotropy, saturation magnetiza¬ 
tion = 8 X 10^ A/m, gyromagnetic ratio 27zy= 1.85 x 10^^ 
rad/T.s and exchange stiffness Aex =13 pJ/m. These values 
of Ms and y were consistent with results from FMR measure¬ 
ments on continuous layers averaged over two reference sam¬ 
ples. The value of y is also close to that determined by Shaw 
et a/®. MNPs were modeled explicitly assuming a ferromag¬ 
netic sphere sitting within the anti-dot with its lower surface 
aligned with the lower surface of the MC. For the MNP, we 
used a = 0.05, A^x and y equal to that in the MC. Ms was 
read off from Fig, [^at the value of F/ext used in the simulation. 
IPython®, Sumatr^SI and matplotlib® were used for analy¬ 
sis, management and visualisation of the simulation data. 

Time domain (ringdown) simulations Time domain or 
Tingdown’ simulations (e.g. Grimsdith et a/.®) were carried 
out using MuMax3® version 3.5.3 with cuboid discretisation 
cells (~ 3.52 X 3.52 x 3.75 nm^) wherein the system’s equi¬ 
librium magnetic configuration, mo(r,F/ext), at a given F/gxt, 
applied in the y-direction, is perturbed with a field pulse in 
the x-direction (0.5 mT sine pulsel^lwith a 300 ps offset and 
30 GHz cut-off frequency). Fourier analysis was then applied 
to the time dependent, spatially averaged x-component of the 
magnetization to extract the characteristic (resonant) frequen¬ 
cies associated with the resultant excited dynamics. Mode vi¬ 
sualisations (insets of Manuscript Fig. 2c) were obtained by 
calculating the spatially resolved Fourier amplitudes for rrix at 
each resonant frequency. 

Eigenmode simulations The second method was an eigen¬ 
mode calculation which uses the finite element micromagnetic 





simulator, FinMag (successor to Nmag^l), to directly deter¬ 
mine the eigenfrequencies and eigenvectors associated with 
a given mo(r,F/ext)- For relaxation, the system was meshed 
using a characteristic internode length of /mesh = 3.5 nm. 
The eigenmodes were determined using an algorithm simi¬ 
lar to that detailed in d’Aquino et al^on a coarsened mesh 
with /mesh = 7 nm (coarsening was needed due to the algo¬ 
rithm’s high memory requirements). The nature of the mode 
(e.g. side or extended) was determined via visual inspection of 
the eigenvectors. 
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